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ABSTRACT 

Based  on  a  new  version  of  Hopcroft  and  Tarjan's  planarity  testing  algo- 
rithm, we  develop  an  O  (mlogn)-nme  algorithm  to  find  a  maximal  planar  sub- 
graph. 
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1.  Introduction 

In  [15],  Wu  defined  the  problem  of  planar  graphs  in  terms  of  the  following  four  subprob- 
lems: 
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Pl .  Decide  whether  a  connected  graph  G  is  planar. 

P2.  Find  a  minimal  set  of  edges  the  removal  of  which  will  render  the  remaining  part  of  G 
planar. 

P3.  Give  a  method  of  embedding  G  in  the  plane  in  case  G  is  planar. 

P4.  Give  a  description  of  the  totality  of  possible  planar  embeddings  of  G  in  the  plane  in 
case  G  is  planar. 

Linear-time  algorithms  for  PI,  P3,  and  P4  have  been  known  for  a  long  time.  The  first 
linear- time  solution  (which  we  call  the  H-T  algorithm)  for  problem  PI  (  the  planarity-testing 
problem)  was  given  by  Hopcroft  and  Tarjan  [7]  in  1974  using  depth-first  search  (DFS)  trees.  A 
P-Q  tree  solution  for  PI  based  on  an  earlier  algorithm  given  by  Lempel,  Even,  and  Cederbaum 
[11]  was  proved  to  have  a  linear-time  implementation  in  1976  partly  by  Even  and  Tarjan  [4]  and 
partly  by  Booth  and  Lueker  [1].  The  P-Q  tree  approach  is  conceptuaUy  simpler,  but  its  implemen- 
tation is  more  complicated  than  that  of  the  H-T  algorithm.  Linear-time  solutions  for  P3  and  P4, 
also  based  on  P-Q  trees,  were  given  by  Chiba  et  al.  [2]  in  1985. 

Wu  [15]  gave  an  algebraic  solution  for  all  four  problems.  He  proved  that  a  graph  is  planar  if 
and  only  if  a  certain  system  of  linear  equations  is  solvable.  In  case  the  graph  is  planar,  an  actual 
embedding  can  be  obtained  by  considering  another  system  of  quadratic  equations.  His  solution  is 
elegant,  but  his  algorithm  takes  O  (m^)  time  on  an  /n-edge  graph. 

Recently,  Jayakumar  et  al.  [9]  studied  problem  P2  ( the  maximal  planar  subgraph  problem). 
For  the  special  case  in  which  a  biconnected  spanning  planar  subgraph  is  given,  their  algorithm 
runs  'mO{n^)  time  and  O  (mn)  space  on  a  graph  with  n  vertices  and  m  edges.  For  more  general 
situations,  their  algorithm  runs  in  0(mn)  time.  Their  algorithm  is  also  based  on  P-Q  trees.  Note 
that  not  every  biconnected  graph  has  a  biconnected  spanning  planar  subgraph  (See  Fig.  1.) 


A  biconnected  graph  that  does  not  have  a 
biconnected  spanning  planar  subgraph 

Fig.  1 


In  this  paper  we  give  an  O  {m\o%n)-\\mt  and  O  (m)-space  solution  to  P2.  For  sparse  graphs 
{i.e.,  graphs  with  m  =  0(jn}^),  where  e  <  1  ),  it  beats  the  algorithm  of  Jayakumar  et  al.  even  in 
the  special  case  when  a  biconnected  spanning  planar  subgraph  is  given.  Indep)endent  of  our  work, 
Di  Battista  and  Tamassia  [3]  have  claimed  an  0(/og/j)-time-per-operation  solution  to  the  prob- 
lem of  maintaining  a  planar  graph  under  edge  additions.  Their  algorithm  also  solves  the  minimal 
planar  subgraph  problem  in  0(m\ogn)  time.  Their  method  is  much  more  complicated  than  ours, 
however,  as  it  is  designed  to  solve  a  more  general  problem.  Recently,  Kobayashi  et  al.  [10]  have 
shown  that  if  a  Hamiltonian  tour  of  the  graph  is  given,  then  P2  can  be  solved  in  liner  lime.  We 
show  that  this  result  can  be  easily  derived  from  our  algorithm  as  a  special  case. 

The  maximal  planar  subgraph  problem  is  closely  related  to  the  planarity-testing  problem.  In 
fact,  a  graph  is  planar  iff  it  is  the  maximal  planar  subgraph  of  itself.  Our  solution  to  the  maximal 
planar  subgraph  problem  is  based  on  the  H-T  algorithm.  But  for  our  purpose,  we  need  to  modify 
the  algorithm.  The  main  difference  is  that  our  version  of  the  algorithm  admits  a  more  general 
ordering  than  the  original  H-T  algorithm  does  in  processing  the  successors  of  each  tree  edge. 
Also,  the  H-T  algorithm  processes  one  path  at  a  time,  while  our  algorithm  processes  one  edge  at  a 
time.  In  this  sense,  our  algorithm  is  a  more  recursive  version  of  the  H-T  algorithm. 

For  the  above  reason,  many  of  our  lemmas  and  theorems  are  similar,  but  not  identical,  to 
those  in  [7].  Instead  of  referring  the  readers  to  [7]  for  the  proofs,  we  find  it  more  convenient  and 
accurate  to  supply  all  main  proofs  in  this  paper. 

The  rest  of  this  paper  is  organized  as  follows.  Section  2  gives  preliminary  definitions.  Sec- 
tion 3  is  a  new  version  of  the  H-T  planarity  testing  algorithm,  which  leads  to  our  maximal  planar 
subgraph  algorithm  in  Section  4.  Section  5  is  a  summary. 

2.  Preliminaries 

Consider  an  undirected  graph  Gq  =  (Yq,  Eq)  with  edge  set  Eq  and  vertex  set  Vq.  Let  /j  = 
I  Vq  I  and  m  =  \Eq\.  We  can  draw  a  picture  Go'  of  Gq  in  the  plane  as  follows:  for  each  vertex 
V  e  Vq,  we  draw  a  distinct  point  v';  for  each  edge  («,  v)  €  Eq,  we  draw  an  simple  arc  connecting 
the  two  points  u  and  v'.  We  call  this  arc  an  embedding  of  the  edge  (u,  v).  For  brevity,  we  will 
sometimes  identify  graphs  with  their  pictures  thus  drawn  on  the  plane.  If  no  arcs  of  Gq'  cross 
each  other,  then  we  call  Gq'  a  planar  embedding,  or  simply  an  embedding,  of  Gq.  If  Gq  has  an 
embedding,  then  we  say  that  Gq  is  planar. 

The  following  facts  are  important  to  our  discussion: 

Observation  1,  Let  C  be  a  simple  closed  curve  in  the  plane  as  in  Fig.  2;  let  a  be  a  point 
inside  C  and  ^  be  a  point  outside  C.  Then  any  curve  that  joins  a  and  b  crosses  C. 
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Fig.  2 


Fig.  3 


Observation  2.  Let  G  i  be  the  undirected  graph  represented  by  Fig.  3,  in  which  P  is  a  path 
joining  the  two  vertices  a  and  b  on  cycle  C.  Then  in  any  embedding  of  G  i ,  all  the  edges  of  path 
P  are  on  the  same  side  of  the  cycle  C  (either  inside  or  outside). 


Fig.  4 


Fig.  5 


Observation  3.  Let  G^  be  the  undirected  graph  represented  by  Fig.  4,  in  which  a\,a2,b\ 
and  ^7  2  are  four  distinct  vertices  that  appear  in  order  on  C.  Then  in  any  embedding  of  G 2,  the  two 
paths  P I  and  P^  are  on  opposite  sides  of  the  cycle  C. 

Observation  4.  Let  G3  be  the  undirected  graph  represented  by  Fig.  5,  in  which  a,  c\,  C2 
and  b  are  vertices  that  appear  in  order  on  C,  and  C\  and  C2  may  be  the  same.  Then  in  any 
embedding  of  G 3,  the  two  subgraphs  P]  and  P2  are  on  opposite  sides  of  the  cycle  C. 

All  four  observations  above  are  intuitively  obvious  and  follow  from  the  Jordan  Curve 
Theorem  [6, 14]. 

A  depth- first-search  {abbr.  DFS)  [7]  will  convert  the  undirected  graph  Go  =  (Vq.  £0)  into  a 
directed  graph  G  =  (V,  T,  B),  where  V  is  the  set  of  DFS  numbers  of  vertices  mVo,T  is  the  set  of 
tree  edges,  and  B  is  the  set  of  back  edges.  Each  edge  of  Gq  is  converted  into  either  a  tree  edge  or 
a  back  edge.  All  the  tree  edges  form  a  DFS  forest.  If  [a,  i>  ]  is  a  tree  edge,  then  a  <  b.  If  [a,  b  ]  is 
a  back  edge,  then  b  <  a,  and  there  is  a  tree  path  in  T  from  ft  to  a.  In  either  case,  a  is  called  the  tail 
of  [a,  b  ],  and  b  is  called  the  head  of  [a,  b  ].  The  union  of  7  and  B  will  be  denoted  by  E. 

For  notational  convenience,  we  will  frequently  identify  undirected  graphs  with  their  DFS 
representations.  Since  we  are  interested  only  in  graphs  with  no  isolated  vertices,  we  will  represent 


graphs  with  their  edge  sets. 

We  define  successors  for  both  vertices  and  edges.  If  [a,  ft  ]  is  a  tree  edge,  then  6  is  a  succes- 
sor of  a.  If  [a,  ft  ]  is  a  tree  edge  and  [ft,  c  ]  is  any  edge,  then  [ft,  c  ]'\s  a  successor  of  [a,  ft].  Back 
edges  have  no  successors.  We  also  define  descendants  and  ancestors  for  both  vertices  and  edges. 
A  descendant  of  vertex  (resp.  edge)  x  is  defined  recursively  as  either  x  itself  or  a  successor  of  a 
descendant  of  x.  If  )>  is  a  descendant  of  x,  then  x  is  an  ancestor  ofy. 

Let  e  =  [a  ,  ft  ]  e  £.  Let  K  be  the  set  of  vertices  y  such  that  for  some  x,  [x,  >> )  is  a  back  edge 
and  also  a  descendant  of  e.  If  Y  is  not  empty,  we  define  lowi(e)  to  be  the  smallest  integer  in  Y, 
and  low2(e)  to  be  the  second  smallest  integer  in  K  u  [n+\ }.  Otherwise,  we  define  lowi{e)  = 
low2{e)  =  n+1.  The  two  mappings  low^  and  lowj  can  be  computed  in  0(m)  time  during  the 
depth- first- search  on  Go  [7].  If  a  is  not  the  root  of  a  DFS  tree,  and  lowi{e)  >  a,  then  a  is  an  arti- 
culation point  of  G  [12]. 

If  e  =  [a,  ft  ]  is  any  edge  in  £,  then  we  define  the  function  (j)  on  £  as  follows. 

\  2  low  lie)  iflow2(e)>a 

[2  low  lie)  +  1     otherwise 

We  arrange  the  successors  of  each  tree  edge  in  increasing  order  on  their  $  values.  This  ord- 
ering can  be  computed  in  0(m)  time  using  a  bucket  sort  [7].  If  ei,  ...,  e^  are  the  successors  of  e 
ordered  this  way,  we  will  call  e,  the  Jth  successor  of  e  for  /  =  1  ..  it. 

As  in  [7],  for  e  =  [a,  ft  ],  we  define  S  (e),  the  segment  of  e,  to  be  the  subgraph  of  G  that  con- 
sists of  all  the  descendants  of  e.  We  use  ATT(e)  to  denote  the  set  of  back  edges  [c,  d]  m  5(e) 
such  that  <i  is  an  ancestor  of  a,  including  a  itself.  Each  back  edge  in  ATT{e)  is  called  an  attach- 
ment of  e. 

For  any  edge  e  =  [a,  ft],  we  define  cycle (e)  as  follows:  if  e  is  a  back  edge,  then  cycle (e)  = 
[e)  u  {e':  e'  belongs  to  the  tree  path  from  ft  to  a};  if  e  is  a  tree  edge  and  low^ie)  >  a,  then 
cycle ie)  =  {);  otherwise,  cycle (e)  =  cycle (ei),  where  e\  is  the  first  successor  of  e.  We  use 
sub(e)  to  denote  the  subgraph  5(e)  u  cycle  (e).  It  is  easy  to  see  that  if  cycle  (e)  is  not  empty,  then 
the  vertex  lowi(e)  is  always  on  cycle (e).  Also,  if  low lie)  >  a,  then  subie)  =  5(e);  if  low^ie)  < 
a,  then  sub (e)  -  5  (e)  =  {e':  e  belongs  to  the  tree  path  from  low  j  (e)  to  a} . 

Fig.  6  illustrates  some  of  these  definitions,  where  lowi(e)  =  1;  low^ie)  =  2;  cycle (e)  =  { 
[1,  2],  [2,  3],  [3,  4],  [4,  5],  [5,  6],  [6,  7],  [7,  8],  [8,  1]  );  5(e)  contains  all  the  edges  in  the  graph 
except  [1,  2],  [2,  3],  [3.  4];  5wft(e)  is  the  whole  graph;  ATT(e)  =  {[8,  1],  [9,  3],  [12,  1],  [14.  2], 
[13,4]. 
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Fig.  6 


3.  Planarity  testing 

As  explained  in  [7],  a  graph  is  planar  if  and  only  if  each  of  its  biconnected  components  is 
planar.  Also,  a  graph  of  one  edge  is  always  planar.  Thus,  we  need  only  consider  how  to  test  the 
planarity  of  biconnected  graphs  with  more  than  one  edge.  Let  G  =  (V,  7,  fi)  be  a  DFS  represen- 
tation of  such  a  graph.  Then  T  forms  a  single  tree  with  only  one  tree  edge  leaving  the  root.  Call 
this  tree  edge  eo-  Since  sub(,eo)  is  the  whole  graph,  we  can  determine  the  planarity  of  G  with  a 
procedure  that  can  determine  the  planarity  of  sub  (e)  for  all  e  in  G. 

We  say  that  an  edge  e  is  planar  if  sub  (e)  is  planar.  To  determine  the  planarity  of  an  edge  e, 
we  consider  two  cases.  If  e  is  a  back  edge,  then  sub  (e)  =  cycle  (e),  which  is  always  planar.  Oth- 
erwise, e  is  a  tree  edge  having  at  least  one  successor.  In  this  case  we  first  determine  the  planarity 
of  each  of  its  successors.  If  all  these  successors  are  planar,  then  we  determine  the  planarity  of  e 
based  on  the  structure  of  its  attachments.  The  details  follow. 


3.1.  Structure  of  attachments 

The  planarity  of  an  edge  e  =  [a,  b]  directly  depends  on  the  suaicture  of  its  attachments. 
Since  we  assume  that  G  is  a  biconnected  graph  with  more  than  one  edge,  then  low^ie)  <  a,  and 
both  ATTie)  and  cycle  (e)  are  not  empty.  If  e  is  planar,  then  we  can  partition  the  edges  of 
ATTie)  into  blocks  as  follows.  We  put  two  back  edges  of  ATTie)  in  the  same  block  if  they  are  on 
the  same  side  of  cycle{e)  in  every  embedding  of  sub (e).  Two  blocks  Bi  and  B2  of  y47T(e) 
interlace  if  they  are  on  opposite  sides  of  cycle (e)  in  every  embedding  of  sub{e).  Each  block  B, 
of  ATT{e)  can  interlace  at  most  one  other  block,  since  two  attachments  of  e  that  cannot  be 


embedded  on  the  same  side  of  cycle  (e)  as  B,  must  be  in  the  same  block. 

The  back  edge  on  cycle  (e)  is  the  only  attachment  of  e  that  will  not  be  embedded  on  either 
side  of  cycle  (e).  By  convention,  this  back  edge  forms  a  block  by  itself,  called  the  singular  block 
of  c,  which  does  not  interlace  other  blocks  oi  ATT(e). 

In  Fig.  6.  ATT(e)  consists  of  four  blocks:  fi,  =  ([8,  1]},  B2  =  {[12.  1],  [14.  2J).  B3  = 
{[9,3]},  and  B  4  =  { [  1 3,  4] } .  B  i  is  singular.  B  2  and  B  3  are  interlacing. 

If  e'  =  [u,  V  ]  is  an  attachment  of  e,  then  low  1  (e)  <  v  <  a.  If  low  i(e)  <  v  <  a,  then  we  say 
that  e'  is  normal.  Otherwise  we  say  that  e'  is  special.  A  block  of  attachments  of  e  is  normal  if  it 
contains  some  normal  attachment  of  e.  Otherwise  we  say  that  it  is  special.  In  Fig.  6.  B 1  and  B4 
are  special,  and  other  blocks  are  all  normal.  We  say  that  sub(e)  is  strongly  planar  w.r.t.  e  if  e  is 
planar  and  if  all  the  normal  blocks  of  ATT (e)  can  be  embedded  on  the  same  side  of  cycle  (e).  If 
sub(e)  is  strongly  planar  (w.r.t.  e),  then  we  say  that  e  is  strongly  planar.  We  have 

Lemma  l.  Let  e,  be  the  ith  successor  of  e,  where  i  >  1 .  Then  e,  is  strongly  planar  iff  the 
subgraph  S  (e,)  u  cycle  (e)  is  planar. 

Proof 

=>  If  e,  is  strongly  planar,  then  there  is  an  embedding  of  sub(ei)  such  that  all  its  normal 
blocks  are  on  the  same  side  of  cycle  (ei).  Thus  we  can  add  cycle  (e)  to  the  other  side  of  cycle  (ei) 
to  get  an  embedding  of  S  (e,)  u  cycle (e). 

<=  If  5(e,)  u  cycle  (e)  is  planar,  then  in  any  embedding  of  5(e,)  u  cycle (e),  all  the  normal 
blocks  of  y47T(e,)  must  be  on  the  same  side  of  cycle  (e,).  ■ 

Note  that  in  an  embedding  of  S  (e,)  u  cycle  (e).  the  special  blocks  of  e,  do  not  have  to  be  on 
the  same  side  of  cycle  (e,).  (See  Fig.  7.) 


cycle  ifii) 


cycle  {e) 


The  two  special  attachments  ct  and  ct'  of  e,  can  be  on  different 
sides  of  cycle  (e,).  although  they  are  on  the  same  side  of  cycle  (e). 

Fig.  7 
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We  will  represent  a  block  of  back  edges  H  =  [[by,  ay],  [b2,  02], ....  [b„  a,]}  by  a  list  K  = 

[fl],  oj flj.  where  ai  <  ^2  -  •••  -  ^i-  Frequently,  we  will  identify  blocks  with  their  list 

representations.  r>ef\ne  first (H)  =  first (K)  =  ay,  and  last {H)  =  last (K)  =  a,.  If  ^T  is  empty,  we 
define  first (H)  =  first (K)  =  n+l,  and  last{H)  =  last{K)  =  0.  We  can  further  organize  the  blocks 
of  ATT (e)  as  follows:  if  two  blocks  X  and  Y  interlace,  we  put  them  into  a  pair  [X,  Y],  assuming 
last  (X)>  last  [Y)\  if  a  nonempty  block  X  does  not  interlace  any  other  block,  we  form  a  pair 
[X,  [  ]].  Let  [Xx,Yx]  and  [X2,  K2]  be  two  pairs  of  interlacing  blocks.  We  say 
[Xx,Yx]<[X2,Y2[\iUast(X\)^min{first{X2),first{Y2)). 

We  say  a  list  of  interlacing  pairs  [q\,...,qs]  is  well-order  ed  xi  q\  <  ■  ■  •  <qs.  Empty  lists  or 
lists  of  one  pair  are  well-ordered  by  convention.  We  will  see  that  all  the  interlacing  pairs  of 
i47T(e)  can  be  organized  into  a  well-ordered  list  [pi,  ...,  p,].  We  call  this  list  a«(e). 

In  Fig.  6,  a«(e)  =  [pi.  P2.P3].  wherep,  =  [[1],  [  ]  ].p2  =  [[3],  [1.  2]],  andpj  =  [[4],  [  ]]. 

3.2.  Computing  on  (e) 

Now  we  are  ready  to  compute  att{e).  The  planarity  of  e  will  be  decided  at  the  same  time. 

Consider  any  edge  e  =  [a,  b[.  If  e  is  a  back  edge,  then  its  only  attachment  is  e  itself.  There- 
fore att(e)  =  [[[b[,  [  ]]].  Otherwise,  let  e\,  ...,  e^  be  the  successors  of  e  in  increasing  order  by 
their  (()  values.  We  first  recursively  compute  att(ei)  for  each  successor  e,  of  e.  Then  we  compute 
att(e)  in  four  steps: 

Algorithm  1 

Step  1  For  i  =  1  ..  k,  delete  all  occurrences  of  b  appearing  in  blocks  within  att{ei).  Because 
these  occurrences  appear  together  at  the  end  of  the  blocks  that  are  contained  in  the  last  pairs  of 
att(ei)  only,  a  simple  list  traversal  suffices  to  delete  all  these  occurrences  in  time  0(k  + 
number  of  deletions).  After  this,  initialize  att  (e)  to  be  att  (e  1 ). 

Step  2.  For  i  -2  ..  k,  merge  all  the  blocks  of  att{ei)  into  one  intermediate  block  B,.  See 
Fig.  8. 

According  to  Lemma  1,  this  step  can  only  be  done  for  a  given  value  of  /  if  the  normal 
blocks  of  att  (e,)  do  not  interlace.  ( If  a  pair  of  normal  blocks  of  att  (e,)  interlace,  the  graph  is  not 
planar,  and  the  compulation  fails.)  To  merge  the  blocks  for  a  given  value  of  i,  we  traverse  the  list 
of  pairs  attie^),  concatenating  blocks  to  form  B,.  Initially,  B,  is  empty.  To  process  pair  [X,  Y],  if 
X  and  Y  are  both  normal,  the  computation  fails,  since  the  graph  is  not  planar.  Otherwise,  we  con- 
catenate X  and  Y  onto  the  end  of  B,  in  order  and  continue.  In  case  that  y  is  a  special  block,  we 
know  that  Y  has  only  attachments  ending  in  lowiiei),  since  b  was  deleted  from  all  blocks  of 
attiCi)  in  Step  1.  Thus  the  correct  ordering  of  attachments  is  maintained  by  this  process.  This 


cycle  (e) 


cycle  (e) 


Fig.  8 


step  takes  0(1+  number  of  blocks  in  an  (e,)),  resulting  in  one  block  for  each  /. 

Step  3.    Merge  blocks  in  att{e).  See  Fig.  9. 

By  Observation  3,  all  blocks  D  in  au{,e)  with  last(D)  >  low  1(62)  must  be  merged  into  one 
block  B I .  (If  any  two  of  these  blocks  interlace,  the  graph  is  not  planar,  and  the  computation 
fails.)  This  is  achieved  by  merging  from  the  high  end  of  att(e).  The  time  is 
0(1  +  reduction  in  number  of  blocks).  This  step  turns  att(e)  into  a  list  of  pairs  pi  <  ■  •■  <Ph 
with  only  pf,  possibly  having  a  block  D  with  last(D)  >  low^ie2).  Note  that  low^iej)  is  the 
lowest  among  the  vertices  low^iej) low\(ek). 


lowx(e2) 


J  ^^  cycle {e) 


Fig.  9 


Step  4.    For  /  =  2, ...,  k,  add  blocks  B,  into  att{e). 

To  process  B,,  consider  the  highest  pair  P:\X,  Y  ]  of  att(,e).  Consider  three  subcases: 
i.  If  B,  cannot  be  embedded  on  either  side  of  cycle  (e),  then  the  computation  of  an  (e)  fails. 
ii.  If  B,  interlaces  X  only,  then  merge  B,  into  Y  by  concatenating  their  ordered  list  represen- 
tations. Next,  switch  X  and  y  if  last  (X)  <  last  (Y). 


10- 


cycle  (e) 


cycle  (e) 


cycle  (e) 


Bi  cannot  be  embedded  in 
either  side  of  cycle  (e) 


S,  interlaces  X  only 
Fig.  10 


Bi  interlaces 
neither  X  nor  Y 


iii.  If  B,  interlaces  neither  X  nor  Y,  then  add  [S,,[  ]]  to  the  high  end  ofattie);  P  :=  [B,,[  ]]. 

By  the  following  lemma,  testing  whether  B,  interlaces  X  or  y  takes  0(1)  time.  Also  by  that 
lemma,  it  is  not  possible  that  B,  interlaces  Y  only,  since  last(X)  >  last{Y)  (  see  Fig.  10 ).      ■ 

Lemma  2.  In  Step  4,  B,  and  D  can  be  embedded  on  the  same  side  of  cycle  (e)  ijf  low  i  (e,)  > 
last  (D).  where  D  =  X  or  D  =  Y. 

Proof 

=>  Assume  lowyie,)  <  last{D).  Then  there  must  he  a  path  P ^  in  S(«,)  from  b  to  /oH'i(e,) 
containing  some  back  edge  in  B, ,  and  another  path  7*2  in  5  (e  i )  u  ...  u  S  (e,_i )  from  a  vertex  w  on 
cycle (e)  to  last{D)  containing  some  back  edge  in  D  but  no  edge  on  cycle {e).  We  consider  two 
cases  (  see  Fig.  1 1).  If  w  >  b,  then  by  Observation  3,  P  \  and  P2  cannot  be  embedded  on  the  same 
side  of  cycle  {e).  If  w  =  i>,  then  the  first  edge  on  P  2  's  Cj  for  some  1  <  y  <  i,  which  implies  /'  >  2 
and  (!)(ej)<())(e,).  Consequently,  low^{ej)  <  low^iei)  <  last{D)  <  b,  which  implies  that 
lowiicj)  <  b.  If  low^iCj)  <  low]{ei),  then  there  must  be  an  undirected  simple  path  P^  between 
last{D)  and  lowx{ej)  containing  back  edges  in  D  but  no  edges  on  cycle{e).  By  observation  3 
again,  P^  and  Pt,  cannot  be  embedded  on  the  same  side  oi cycleie).  \i  low^icj)  =  lowi(e,),  then 
/oH'2(e,)  <  b  (recall  lowiicj)  <  low](e,)).  By  observation  4,  S(e,)  and  S(ej)  cannot  be  embedded 
on  the  same  side  of  cycle  (e).  All  of  the  above  cases  imply  that  B,  and  D  cannot  be  embedded  on 
the  same  side  of  cycle  (e). 

<=  See  the  proof  of  Lemma  4  in  the  next  section.      ■ 


THEOREM  1. 
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■ycleie)         /owjCe,) 


w  =6  and 
low  1  (e,)  >  low  1  (Cj) 

Fig.  11 


"cycle  (e) 


w  =  b  and 


1.  Algorithm  1  computes  att{e)  successfully  iff  e  is  planar. 

2.  If  e  is  planar,  then  Algorithm  1  computes  attie)  correctly. 
Proof    See  the  next  section.  ■ 


33.  Correctness 

In  the  following  proofs,  unless  slated  otherwise,  we  will  use  att(e)  to  mean  the  list  att{e) 
computed  by  Algorithm  1.  But  we  will  prove  that  this  att{e)  correctly  implements  the  att(e) 
defined  in  Section  3.1. 

During  the  presentation  of  Algorithm  ] ,  we  explained  that  two  nonempty  blocks  form  a  pair 
within  att(e)  only  if  they  cannot  be  embedded  on  the  same  side  of  cycle  (e),  and  the  computation 
of  att(e)  fails  only  when  e  is  not  planar.  Also  we  can  see  that  the  singular  block  of  e  is  not 
merged  with  any  other  block.  To  prove  Theorem  1,  we  still  have  to  show  that  the  following 
assertions  are  true: 

(1)  if  computation  of  att{e)  succeeds,  then  e  is  planar, 

(2)  if  any  two  nonempty  nonsingular  blocks  of  att{e)  do  not  form  a  pair,  then  these  blocks 
can  be  embedded  on  the  same  side  as  well  as  on  different  sides  of  cycle  (e); 

(3)  attie)  is  well-ordered. 
We  first  prove  (3),  i.e.. 

Lemma  3.  The  list  of  pairs  att(e)  computed  by  Algorithm  1  is  well-ordered. 

Proof  We  prove  this  lemma  by  induction  on  the  number  of  descendants  of  e.  If  e  has  no 
successor,  then  e  is  a  back  edge,  and  the  lemma  is  trivially  true.  Now  assume  that  e  is  a  tree  edge 
with  successors  e\,  ...,  e*  in  increasing  order  by  $  value,  and  that  att{e\), ....  attie^)  are  all  well- 
ordered.  After  Steps  1  and  2  are  executed,  off  (e  i) att{ek)  are  still  well-ordered.  Thus,  attie) 
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is  well-ordered  when  it  is  initialized  to  att(ei).  In  Step  3,  only  blocks  in  the  highest  pairs  of 
att(e)  are  merged,  and  therefore  att(e)  is  still  well-ordered  after  the  merge.  Then  consider  the 
moment  in  Step  4  just  before  B,  is  added  to  att{e).  Assume  arr(e)  is  well-ordered  at  this  moment. 
Let  P.[X,  y]  be  the  last  pair  of  att(e).  We  need  only  consider  the  two  cases  in  which  the  com- 
putation does  not  fail. 

1.  last(Y)  <  lowiiBi)  <  last(X).  Then  fi,  is  merged  with  Y.  If  P  is  the  only  pair  in  att(e), 
then  att(e)  is  well-ordered  by  definition.  Otherwise,  let  Q:  [X\,  Ki]  be  the  pair  next  to  P.  Then 
we  have  Q  <P  before  merge.  We  need  only  to  show  that  this  is  still  true  after  the  merge.  If  /  =  2, 
then  Step  3  guarantees  that  first(B2)  =  lowi(e2)  ^  max(last(Xi),  last(Y\)).  Hi  >  2,  then  fi,_i  is 
contained  in  either  X  or  Y.  Since  first (Bi)  =  lowi(ei)  >  lowi(ei-\)  =  first {Bi-\)  > 
min{first{X),  first (Y)),  then  merging  B,  into  Y  does  not  change  the  value  of 
mxnifirstiX),  first  (Y)).  Thus,  after  merging  Y  and  B,,  we  still  have  Q<P. 

2.  lowiiBi)  >  last(X).  Then  [B,,  [  ]]  becomes  the  last  pair  of  att{e).  Since  last(X)  < 
low  \{ei)=  first  {Bi)  in  this  case,  we  have  P  <  [B,,  [  ]]. 

Thus,  att(e)  is  still  well-ordered  after  each  B,  is  added,  i  =2  ..  k.  Therefore  att(e)  is  well- 
ordered  after  Step  4.  ■ 

Next  we  prove  the  j/part  of  Lemma  2:  In  Step  4  of  Algorithm  1,  if  low  ^(ei)  >  last(D),  then 
B,  and  D  can  be  embedded  on  the  same  side  of  cycle (e),  where  D  =Xot  D  =  Y. 

Proof  Consider  an  embedding  of  sub{e)  before  B,  is  added.  Assume  without  loss  of  gen- 
erality that  X  is  embedded  on  the  left  hand  side  of  cycle  {e)  in  this  embedding.  Let  h  =  last{X). 
Then  h  =  max(/asr(Z):  Z  is  a  block  in  ««,),  and  there  is  a  face  F  on  the  left  of  cycle (e)  in  the 
current  embedding  of  sub  (e)  such  that  the  tree  path  from  h  to  bis  on  the  boundary  of  F.  See  Fig. 
12. 


cycle  (e) 
Fig.  12 

Thus,  if  low liBi)  >  last(X),  then  B,  can  be  embedded  in  F.  Similarly,  let  h'  =  ma\{last(Uy.  U  is 
block  in  atti  embedded  on  the  right  hand  side  of  cycle  (e)}.  Then  there  is  a  face  F'  on  the  right  of 
cycle (e)  in  the  current  embedding  of  sub{e)  such  that  the  tree  path  from  h'  to  b  is  on  the 
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boundary  of  F.  According  to  the  proof  of  Lemma  3,  if  low^ifii)  >  last{Y),  then  /owiCB,)  >  h' 
also.  Therefore  B,  can  be  embedded  in  F'  in  this  case.  ■ 

Then  we  prove  the  following  lemma  that  implies  the  assertions  (1)  and  (2).  We  say  that  a 
set  W  of  blocks  of  att{,e)  is  consistent  w.r.t.  e  if  for  aU  X,  K  e  W,  neither  [X,  Y]  nor  [Y,  X]  is  in 
au{e). 

LEMMA  4..  If  Algorithm  1  does  not  fail,  and  D  j  and  D2  are  two  disjoint  consistent  sets  of 
nonsingular  blocks  from  att(e),  then  there  is  an  embedding  of  sub  (e)  such  that  blocks  of  D 1  are 
on  one  side  of  cycle  (e)  and  blocks  ofDi  are  on  the  other  side  of  cycle  (e). 

Proof    The  lemma  is  trivially  true  if  e  has  no  successors.  If  e  has  successors,  let  e  1 e* 

be  the  list  of  successors  of  e  in  increasing  order  by  their  ({)  values.  Assume  that  the  lemma  holds 
for  each  of  these  successors.  We  want  to  construct  an  embedding  of  sub  {e)  such  that  D  j  and  D^ 
are  embedded  on  different  sides  of  cycle  (e). 

If  H'  is  a  set  of  blocks,  then  a  W-attachment  is  an  attachment  contained  in  some  block  of  W. 
Forj=  1.2.  letD/=  [X:[X,  y]or[y.  X]  is  a  pair  inaff(e)  and  y  e  Dj}.  Let//]  =  (  X:  ([X,  Y] 
is  a  pair  in  att(e))  and  (  X  and  y  are  not  in  D ,  u  D2)  ).  and  let  Hj  =  {  Y:  ([X,  Y]  is  a  pair  in 
attie))  and  (X  and  y  are  not  in  D,  ^  Dj)  )•  Let  Ci  =0,  UD2'  uHi,  and  C2=D2  uD/u 
H2.  For;  =  1,2,  let  Kj  =  {X:  X  is  a  block  in  5a/7 (e  1 )  containing  some  Cj -attachment}.  Then  Ki 
and  K2  are  two  disjoint  consistent  subsets  of  blocks  of  sub  (e  \ ). 

Initially,  we  construct  an  embedding  of  subieO  such  that  A^i  and  K2  are  on  different  sides 
of  cycle  (e  1 ).  As  a  result,  those  C 1  -attachments  and  C  2  -attachments  contained  in  sub  (e  1 )  are  on 
different  sides  of  cycle {e)  (which  is  cycle {e\)).    This  embedding  exists  by  the  induction 

hypothesis.  Take  this  embedding  to  be  the  initial  embedding  of  att{e).  Then  for  i  =  2 k,  wc 

add  sub  (e,)  to  this  embedding  one  by  one  as  follows. 

Since  the  normal  blocks  of  arr  (e,)  do  not  interlace,  we  can.  by  induction,  find  an  embed- 
ding of  subiei)  such  that  all  of  its  normal  blocks  are  embedded  on  the  same  side  of  cycle  (ei).  We 
call  this  embedding  £,.  and  its  mirror  image  £,'.  Let  fi,,  /*  and  [X,  y  ]  be  the  same  as  in  Step  4  of 
algorithm  1.  Assume  without  loss  of  generality  that  X  is  embedded  on  the  left  hand  side  of 
cycle  (e).  Consider  two  cases: 

Case  1:  last(,Y)  <  /owjCe,)  <  last(X).  Then  B,  is  merged  with  Y.  According  to  Lemma  2, 
one  of  £,  or  £,'  can  be  embedded  on  the  right  of  cycle  (e).  If  B,  contains  any  C 1  -attachment,  then 
X  contains  some  C 2 -attachment;  if  B,  contains  any  C 2 -attachment,  then  X  contains  some  C\- 
attachment.  In  any  case.  C 1 -attachments  and  C  2 -attachments  are  still  on  different  sides  of 
cycle  (e)  after  B,  is  embedded. 

Case  2:  tow  ,(e,)  >  last(X).  Again  by  to  Lemma  2.  one  of  £,  or  £,'.  say  E,.  can  be  embed- 
ded on  the  left  of  cycle  {e),  and  the  other.  £,'.  can  be  embedded  on  the  right  of  cycle  (e).  One  of 
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these  two  choices  will  result  in  an  embedding  such  that  C  i -attachments  are  on  one  side  of 
cycle  (e),  and  C 2 -attachments  are  on  the  other  side. 

When  all  the  B,'s  are  added,  we  get  an  embedding  of  sub  (e)  such  that  D 1  is  on  one  side  of 
cycle  (e),  and  D2  is  on  the  other  side  of  cycle  (e).  TTiis  is  true  because  every  D 1 -attachment  is  a 
C 1 -attachment,  and  every  D2-attachment  is  a  C 2 -attachment.  ■ 

This  completes  the  proof  of  Theorem  1,  and  establishes  that  the  list  att(e)  computed  by 
Algorithm  1  has  the  properties  discussed  at  the  end  of  Section  3. 1. 

Let  e,  e,,  fi,,  arr,,  X,  Y,  and  h'  be  the  same  as  in  the  above  proofs.  Let  L  =  [Z:  [Z,  f/]  is  a 
pair  in  att,],  R  =  {U:  [Z,  U]  \s  a.  pair  in  an,),  and  hj  =  {last{U):  f/  is  a  block  in  R].  Then  it  is 
easy  to  see  that  h'  >  /12.  According  to  Lemma  4,  there  exists  a  embedding  of  sub{e)  (before 
adding  B,)  such  that  L  is  embedded  in  one  side  of  cycle  (e)  and  R  is  embedded  on  the  other  side. 
According  to  the  proof  of  Lemma  2,  low ,  (e,)  >  last  (Y)  iff  low  1  (e,)  >  /i  2 .  Therefore  we  also  have 

Corollary  l .  l .  B,  cannot  be  embedded  in  either  side  of  cycle  (e)  iff  low  j  (e, )  <  /i  2  •    ■ 
Corollary  1 . 1  gives  a  test  whether  B,  can  be  added  to  atti  without  referring  to  the  top  pair  of 

atti.  This  is  useful  in  our  maximal  planar  subgraph  algorithm,  where  we  need  to  test  whether  B  2 

can  be  added  to  att(e)  before  Step  3  is  performed. 

3.4.  Data  structure  and  running  time 

As  suggested  in  [7],  we  can  implement  blocks  as  linked  lists.  An  interlacing  pair  of  blocks 
can  be  represented  as  a  record  containing  two  pointers  to  the  two  linked  lists  representing  these 
two  blocks.  Then  att(e)  can  be  represented  as  a  linked  list  of  such  records.  In  this  way,  the  time 
cost  for  Step  1  is  Oik  +  number  of  deletions).  The  cost  for  Step  2,  3  and  4  is 
Oik  +  reduction  in  number  of  blocks).  The  cumulate  expense  of  executing  these  steps  over  the 
whole  graph  is  0(m).  The  initial  DFS  in  which  /ow,  values  are  computed  takes  time  0(m). 
Arranging  the  successors  in  increasing  order  by  (J)  value  for  all  tree  edges  takes  0(m)  time  using  a 
bucket  sort.  Thus  the  whole  algorithm  runs  in  Oim)  time.  It  is  well  known  that  any  (9(m)-time 
algorithm  for  planarity  testing  can  be  implemented  in  0{n)  time  since  m  =  0{n)  for  a  planar 
graph  [7]. 

3^.  A  modification  to  Algorithm  1 

Consider  Step  4  of  Algorithm  1 .  Lemma  2  requires  that  the  successors  of  each  tree  edge  be 
ordered  by  (j)  values.  Maintaining  this  ordering  causes  difficulties  in  solving  the  maximal  planar 
subgraph  problem.  Fortunately,  we  can  modify  Algorithm  1  so  that  it  requires  only  the  low  1  ord- 
ering of  the  successors  of  each  tree  edge. 
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Let  e  =  [a,  fc]  be  a  tree  edge,  and  e  \ e*  be  the  list  of  its  successors  in  increasing  order 

by  low]  values.  Still  define  cycle (e)  =  cycle (ci).  Then  Step  1,  2,  and  3  can  be  perfomed  w.r.t. 
this  ordering  without  any  modification. 

Next  we  want  to  merge  B2,  ....  fl*  into  att(e)  in  that  order.  In  general,  successors  ordered 
by  low  1  values  may  not  be  ordered  by  (J)  values.  Consequently,  there  may  be  some  1  <  i  <k  such 
that  <t>(e,_i )  >  <t>(e,).  But  if  this  happens,  we  know  that  low  1  (e,_i )  =  low  1  (g,)  and  /ovv2(e,)  >  b.  If 
i  =  2,  Lemma  2  still  applies,  and  we  can  merge  B2  into  a«(e)  as  before.  Otherwise,  the  following 
lemma  says  that  we  do  not  have  to  merge  B,  into  att(e)  at  all: 

LEMMA  5.  If  for  some  2  <  i  <  k,  /ow  i(e,_i)  =  lowiiei),  low2(ei)>b,  and  ei  is  planar,  then 
G  is  planar  iJfG-S  (e,)  is  planar. 

Proof  The  only  if  pan  is  trivial,  so  we  just  prove  the  if  pan.  Consider  an  embedding  £,  of 
G,  =G  -5(e,).  Under  the  condition  of  the  lemma,  e,  has  no  normal  attachments.  Since  e,  is 
planar,  then  e,  is  strongly  planar.  Also,  b  and  /owi(e,)  are  the  only  two  vertices  shared  by  5(e,) 
and  G,.  Therefore  S(e,)  can  be  embedded  in  any  face  of  £,  whose  boundary  contains  the  two  ver- 
tices 6  and  lowiici). 

Let  P  be  the  tree  path  cycle (e,-i)  n  cycle {e)  and  let  C  be  the  closed  curve  cycle ie,-i)  u 
cycle  (e)  -  P.  Then  C  contains  edges  from  both  S (e,_i )  and  G,  -  5 (e,_i ).  By  Observation  2,  P  is 
on  one  side  of  C.  Call  this  side  of  C  Sj,  and  the  other  side  S2.  Let  {/  be  the  set  of  faces  in  ^2 
whose  boundaries  contain  edges  firom  5(e,_i)  only,  and  let  W  be  the  set  of  faces  in  Sj  whose 
boundaries  contain  edges  from  G,  -5 (<?,_])  only.  Then  faces  in  U  and  faces  in  W  do  not  share 
common  boundaries.  Thus,  within  5  2  there  must  be  some  face  F  whose  boundary  contains  edges 
from  both  5(e,_i)  and  G,-5(e,_i),  and  therefore  contains  at  least  two  vertices  common  to 
S (e,_i )  and  G,  -  5 (e,_i ).  But  all  the  vertices  common  to  S (e,_i )  and  G,  -  S (e,_i )  are  on  P,  and 
among  them  only  b  and  lowi(e,-\)  are  on  the  boundary  of  52-  Therefore  these  two  vertices  must 
be  on  the  boundary  of  F.  Thus  we  can  embed  S (e,)  in  F  to  get  an  embedding  of  G.  ■ 

Therefore,  under  the  conditions  of  Lemma  5,  in  deciding  the  planarity  of  G,  we  can  ignore 
its  subgraph  5(e,).  Since  the  condition  /ow ](«,_])  =  toH'i(e,)  and  low2(ei)>b  is  implied  by 
low  1  (e,_i )  <  low  1  (e,)  and  <t)(e,_i )  >  <t)(e,),  we  can  modify  Step  4  as  follows: 

Step  4'.   Add  blocks  B2 B*  into  attie)  in  that  order,  assuming  lowi(ei)  <lowi(e2) 

<  •  •  •  <  lowiick).  Initially,  let  j  =  1  and  /  =  2.  To  process  B,,  we  consider  two  cases.  U  j  =  1  or 
<!)(«;)  <  (|)(e,),  we  do  the  same  thing  as  in  Step  4,  and  then  let  J  =  i;  otherwise,  we  do  nothing. 

The  list  att(e)  computed  by  the  modified  algorithm  may  not  contain  all  the  attachments  of 
e.  Some  attachments  may  be  omitted  by  Step  4',  because  their  existence  does  not  affect  the 
planarity  of  the  whole  graph  G. 
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4.  The  maximal  planar  subgraph  problem 

Now  we  consider  the  maximal  planar  subgraph  problem:  find  a  minimal  set  of  edges  whose 
deletion  results  in  a  planar  graph.  The  resulting  graph  is  called  a  maximal  planar  subgraph  of  G. 
We  can  always  find  a  maximal  planar  subgraph  of  G  by  deleting  back  edges  only,  since  all  the 
tree  edges  form  a  forest,  which  is  planar. 

We  will  not  assume  that  the  input  graph  is  biconnected,  since  deletion  of  back  edges  may 
turn  a  biconnected  graph  into  a  graph  with  articulation  points.  But  without  loss  of  generality  we 
can  assume  that  the  input  graph  is  connected.  Thus  the  tree  edges  of  G  form  a  single  tree  with 
root  r.  Let  ri , ....  r^  be  the  tree  edges  leaving  the  root.  If  5  =  1,  then  sub(ti)  is  the  whole  graph  G. 
If  .s  >  1,  then  r  is  the  only  vertex  common  to  subOi),  ••-  sub{ts).  Thus,  to  find  a  maximal  planar 

subgraph  of  G,  we  can  just  find  a  maximal  planar  subgraph  for  each  of  the  subgraphs  sub  (r  j ) 

sub  its),  and  then  simply  put  these  subgraphs  together.  Therefore,  what  we  need  is  a  procedure 
that  can  find  a  maximal  planar  subgraph  of  sub  {e)  for  any  given  edge  e  of  G. 

4.1.  Maximal  l-planar  subgraphs 

We  cannot  build  a  maximal  planar  subgraph  of  sub  {e)  by  simply  putting  together  the  maxi- 
mal planar  subgraphs  of  sub{e\),  ...,  sub{ei,),  and  deleting  those  back  edges  causing  failure  in 
Algorithm  1.  The  reason  is  that  after  these  edges  are  deleted,  it  may  turn  out  that  some  other 

edges,  which  we  deleted  for  making  sub{e\) sub{ek)  planar,  would  not  have  had  to  be  deleted 

at  all.  We  avoid  this  difficult  situation  by  constructing  such  maximal  subgraphs  S\ Sk  of 

sub{e\),  ....  subifik)  that  they  can  be  used  to  construct  a  planar  subgraph  5  of  subie)  without 
further  deletion  of  edges.  Two  measures  are  taken  for  this  purpose.  Firstly,  those  back  edges  in 
sub{e,)  that  can  cause  failure  in  Step  3  or  Step  4  of  Algorithm  1  are  deleted  before  a  maximal 
subgraph  of  S,  is  recursively  computed.  Secondly,  the  information  where  blocks  of  subici)  are 
allowed  to  interlace  is  passed  to  the  recursive  call  that  computes  5,,  so  that  when  the  returned  S, 
is  merged  to  sub{e).  Step  2  of  Algorithm  1  can  also  be  performed  successfully  without  deletion. 
Since  the  planar  subgraph  S  of  sub{e)  computed  by  our  algorithm  may  be  used  to  build  a  larger 
planar  subgraph  of  G  in  the  same  way  as  we  use  S\,  Sj,  ■-,  S^  to  build  S,  we  also  need  to  know 
where  in  S  blocks  are  allowed  to  interlace.  This  approach  leads  naturally  to  the  concept  of 
l-planar  subgraphs,  which  is  a  generalization  of  the  concept  of  strongly  planar  subgraphs. 

Consider  an  edge  e  =  [a,  b]  and  a  vertex  /  on  the  tree  path  from  lowi(e)  to  a.  An  attach- 
ment [u,  v]  of  e  is  l-normal  if  lowx{e)  <  v  <  /.  A  block  of  attachments  is  l-normal  if  it  con- 
tains some  /-normal  attachment.  Let  D  be  the  list  representation  of  a  nonempty  block  of  attach- 
ments. Define  second{D)  to  be  the  second  smallest  element  in  the  set  [x:  j:  e  D)  u  {n+\ ),  and 
define  second {[Y)  =  n+\.  Then  D  is  /-nomial  iff  low-^ie)  <first{D)  <  I  or  second{D)  <  I.  The 
two  mappings  first  and  second  can  be  maintained  during  the  computation  of  attie)  in  0(1)  time 
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for  each  modification  to  att(e). 

We  say  that  the  subgraph  sub(e)  is  l-planar  if  e  is  planar  and  the  /-normal  blocks  oi  att{e) 
do  not  interlace,  (c./.  Fig.  13,  where  (a)  is  /-planar,  but  (b)  is  not)  Edge  e  is  l-planar  \f  sub(e) 
is  /-planar. 

If //  is  a  subgraph  obtained  from  sub{e)  by  deleting  back  edges  only,  then  we  can  define  the 
/-planarity  for  H  (w.r.t.  e)  in  the  same  way  as  we  did  for  sub{e).  We  will  talk  about  /-planar  sub- 
graphs of  sub{e)  in  this  sense.  An  /-planar  subgraph  of  sub{e)  is  maximal  if  it  can  be  obtained 
from  sub(e)  by  deleting  a  minimal  set  of  back  edges. 

Consider  edge  e  =  [a,  b].  According  to  our  definition,  e  is  planar  iff  e  is  /ovi'i(e)-planar, 
and  e  is  strongly  planar  iff  e  is  a-planar.  Therefore,  if  we  can  find  a  maximal  /-planar  subgraph  of 
sub  (e)  for  any  /  with  low  j  (e)  <  /  <  a,  then  we  can  compute  a  maximal  planar  subgraph  of  sub  (e). 


cycle  (e) 


cycle  {e) 


(a) 


(b) 


Fig.  13 


The  following  is  an  outline  of  our  maximal  /-planar  subgraph  algorithm,  where  /  is  a  given 
integer  with  lowi{e)<  I  <a  and  remains  fixed  during  the  processing  of  an  edge.  Let  e  =  [a,  b], 
and  consider  three  cases: 

Case  1 :  e  is  a  back  edge.  Assign  [[[b  ],[]]]  to  att(e),  and  return. 

Case  2:  e  is  a  tree  edge  with  no  successors.  Assign  [  ]  to  att{e),  and  return. 

Case  3:  e  is  a  tree  edge  with  successors  e\ e^,  among  which  e\  has  the  smallest  lowx 

value.  We  construct  a  sequence  G\,  ...,  G^  of  /-planar  subgraphs  of  sub(e)  such  that  Gi  is  a 
maximal  /-planar  subgraph  of  subiC])  and  /owi(ei)  remains  unchanged;  G^  is  a  maximal  /- 
planar  subgraph  of  subie);  and  each  G,,  ]  <  i  <  k,  is  obtained  from  G,_i  by  adding  to  it  a 
strongly  planar  subgraph  5,  of  sub{ei),  where  e,  is  some  successor  of  e  not  contained  in  G,-\. 
During  the  construction,  we  compute  att{e)  using  the  modified  version  of  Algorithm  1.  We 
describe  below  in  rough  terms  how  we  compute  S, : 

1         select  an  edge  e,  with  the  smallest  low  j  value  from  successors  of  e  not  contained  in  G,_i ; 
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2  while  there  exists  a  maximal  strongly  planar  subgraph  of  sub  (e,)  whose  addition  to  G,_i 

destroys  its  /-planarity  do 

3  delete  some  attachments  from  sub(ei)\ 

4  if  the  deletion  changes  the  low  i  value  of  e,  then 

5  select  a  possibly  new  edge  e,  with  the  smallest  low  j  value  from 

successors  of  e  not  contained  in  G,_i ; 

6  end  if; 

7  od; 

8  recursively  construct  a  maximal  strongly  planar  subgraph  of  sub  (e,)  without  changing 
low\(ei)  further.  We  take  this  subgraph  as  5,. 

In  the  procedure  sketched  above,  lines  1,  4,  5,  6,  and  8  guarantee  that  subgraphs  S,  are  gen- 
erated in  increasing  order  by  new  low^  values  of  the  corresponding  successors.  For  each 
1  <  /  <  *,  once  Si  is  computed,  no  edges  will  be  deleted  further  from  it.  There  are  still  two  ques- 
tions remaining  to  be  answered:  how  the  testing  in  line  2  can  be  done  without  constmcting  a  max- 
imal strongly  planar  subgraph  of  sub  (ci),  and  how  the  attachments  are  chosen  so  that  the  deletion 
in  line  3  makes  the  set  of  deleted  edges  minimal.  These  two  questions  are  closely  related  and  will 
be  explained  together  in  the  next  section. 

Remark  In  Algorithm  1,  we  do  not  need  the  concept  of  /-planarity,  since  our  purpose  is  to 
check  the  planarity  of  G.  If  some  interlacing  blocks  of  sub  (e)  are  found  not  being  able  to  fit  in  the 
whole  graph  after  returning  from  several  levels  of  recursive  calls,  we  simply  declare  that  the 
graph  is  not  planar.  But  if  we  want  to  construct  a  maximal  planar  subgraph  of  G,  then  it  is  too 
late  to  delete  edges  efficiently  by  that  time.  Therefore  we  use  the  parameter  /  to  pass  the  informa- 
tion where  blocks  of  sub{e)  are  allowed  to  interlace,  to  the  recursive  calls,  so  that  the  correct 
edges  are  already  deleted  during  the  processing  of  sub  (e). 

The  need  to  generalize  to  /-planarity  arises  in  the  following  way  in  the  algorithm  sketched 
above.  To  compute  a  maximal  planar  subgraph  of  the  input  graph,  the  recursive  calls  that  con- 
struct S,  for  /  =  2, ....  k  must  construct  maximal  strongly  planar  graphs.  Within  one  of  these  recur- 
sive calls,  the  initial  second-level  recursive  call  (to  construct  a  maximal  6-planar  subgraph  of 
sub{ei\),  where  en  is  the  first  successor  of  edge  e,)  and  more  deeply  nested  recursive  calls  of  the 
same  kind  construct  maximal  /-planar  subgraphs  for  general  values  of  /.       ■ 

4J,.  Algorithm  for  deleting  back  edges 

Let  e  =  [a,  b],  and  consider  the  while-loop  in  the  procedure  sketched  above.  If  lowiie,)  > 
b,  then  b  is  the  only  vertex  common  to  subici)  and  G  -sub(e,).  In  this  case,  we  can  apply  the 
maximal  planar  subgraph  to  sub{e,)  separately,  and  do  not  have  to  consider  the  effect  on  the 
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whole  graph.  Next,  we  consider  the  case  when  low\(ei)  <  b.  Assume  that  sub(ei)  is  made 
strongly  planar  by  deleting  some  back  edges.  Suppose  that  the  low^  value  and  the  I0W2  value  of 
e,  are  not  changed  by  these  deletions.  We  want  to  see  whether  the  union  of  subiei)  and  G,_i  is 
/-planar. 

As  in  planarity  testing,  let  5,  be  the  block  of  attachments  obtained  by  merging  att(e,y,  let 
atti  be  the  current  value  of  att{,e)\  let  Bj  be  the  last  block  merged  into  a«,  by  Step  4';  let  /i  1  = 
max{/ajr(Z):  [Z,  i/]  is  a  pair  in  arr,}  and  /12  =  md.\[last{U):  [Z,  i/]  is  a  pair  in  atti).  (Initially, 
we  set  j  =1  and  attie)  =  att(e\)  after  removing  all  the  occurrences  of  b  from  an(ei).)  Finally, 
let  /13  =  max.[last(Zy.  Z  is  an  /-normal  block  of  atti]. 

The  two  variables  h\,  hi  can  be  maintained  in  0(1)  time  per  modification  to  attie)  by 
maintaining  two  lists  L  and  /?  (  as  suggested  in  [7]  ),  where  L  is  the  ordered  list  of  nonempty 
blocks  X  such  that  [X,  K]  is  a  pair  in  att(e),  and  R  is  the  ordered  list  of  nonempty  blocks  Y  such 
that  [X,  y]  is  a  pair  in  att{e).  U  Bi  and  Bn  are  the  highest  blocks  of  L  and  R  respectively,  then 
/i]  =  last{Bi)  and  hi  =  lastiBn).  Lists  L  and  R  also  let  /23  be  maintained  easily.  If  e  is  a  back 
edge,  /13  =  0  from  its  definition.  If  e  is  a  tree  edge,  we  get  the  initial  value  of  /13  from  the  compu- 
tation of  att{e]),  and  modify  it  in  0(1)  time  for  each  modification  of  att(e).  The  details  will  not 
be  discussed  here. 

By  Lemma  5,  in  case  that  e,  is  strongly  planar,  sub(ei)  can  affect  the  planarity  of  G  only  if 
any  of  the  following  conditions  holds; 

a.  7  =  1,  i.e.,  no  block  Bj  has  been  merged  into  attie)  yet, 

b.  low  1  iej)  <  low  1  ie,),  or 

c.  lowiiei)  <  b. 

If  any  of  these  conditions  is  true,  we  consider  two  additional  cases: 

1.  The  union  of  subiei)  and  G,_i  is  not  planar.  By  Corollary  1.1,  this  happens  iff 
lowxici)  <  /l2- 

2.  The  union  of  subie^)  and  G,_i  is  planar,  but  not  /-planar.  Then  6,  is  /-normal,  and  it 
interlaces  an  /-normal  block  of  att,.  We  know  that  B,  is  /-normal  iff /o>vi(ei)  <  low^iei)  <  I  or 
lowiiei)  <  I.  Also,  it  is  easy  to  see  that  B,  interlaces  an  /-normal  block  of  att^  iff  low  ]  (e,)  <  /i  3 . 

This  means,  under  the  conditions  a,  b,  ore,  that  the  union  of  sub  ie,)  and  G,_i  is  not  /-planar 
iff  any  of  the  following  conditions  holds 

i.  low^iei)  <  hi 

ii.  lowiici)  <lowiie,)<  min  {/13,  /} 

iii.  lowiiei)  <  /  and  /ow](e,)  <  /13 
Therefore,  if  any  of  the  condition  i,  ii,  and  iii  holds,  some  back  edge  has  to  be  deleted  to  change 
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either  the  lowi  value  or  the  I0W2  value  of  e,.  Such  testing  and  deletion  can  be  done  even  before 
making  sub(e,)  strongly  planar.  For  this  purpose,  we  combine  the  above  conditions  (a  or  b  or  c) 
and  (i  or  ii  or  iii)  into  two  groups  according  to  whether  they  involve  low-iiei)  or  not: 
Condition  AA. 

(j  =  1  and  low  j  (e,)  <  A 2)  or 

(lowiicj)  <  lowi(ei)  <  h2)  or 

{low\{ej)  <  lowi(ei)  <  min{/j3,/}) 
Condition  BB. 

(/ow2(e,)  <  b  and  /owi(e,)  <  hi)  or 

(low2(ei)  <  band  lowi(e\)  <  lowi(ei)  <mm{h^,  /})or 

(/oH'2(e,)  <  /  and  /ow i(e,)  <h^) 

It  can  be  checked  that  the  condition  ((a  or  b  or  c)  and  (i  or  ii  or  iii))  is  equivalent  to  the  con- 
dition (AA  or  BB). 

If  Condition  AA  is  true,  we  can  make  it  false  only  by  changing  the  value  /owjCe,).  In  this 
case,  we  delete  all  the  back  edges  of  sub  (Ci)  entering  the  vertex  towiCe,).  After  the  deletion,  we 
choose  a  pxjssibly  new  e,  with  the  smallest  low  j  value. 


low  I  (Cj 


cycle  (e) 


The  edge  e,  satisfies  condition  BB.  If  we  choose  to  delete  cf,  then  d"  wiU  also 
be  deleted  later  because  of  Condition  AA,  and  the  resulting  graph  will  not  be  maximal. 

Fig.  14 


If  Condition  AA  is  false,  then  we  test  Condition  BB.  If  the  result  is  true,  we  know  that 
/owi(e,)  =  /oH'i(e^).  This  is  because  /ovvjCe,)  <  lowiici),  which  means  that  BB  implies  that 
lowiici)  <  hj  or  loWi(ei)  <  mm[h^,  /),  from  which  it  follows  that  lowi(ej)  =  /owjCe,);  otherwise 
AA  would  be  true.  (We  have  low^iej)  <  /owi(e,)  by  the  ordering  of  the  successors  of  e.)  To 
make  Condition  BB  false,  we  can  change  the  value  of  either  low  j  (e,)  or  lowiiei).  If  we  choose  to 
change  /ow2(e,)  consistently,  then  at  least  one  of  the  back  edges  [u,  v]  of  subici)  with 
V  =  low  1  (e,)  will  survive.  But  if  we  choose  to  change  low  1  (e,),  it  may  happen  that  all  the  attach- 
ments in  ATT(ei)  n  A7T(e)  are  deleted  eventually  and  that  the  resulting  graph  is  not  maximal. 


-21  - 

Therefore,  in  this  case  we  choose  to  delete  all  the  back  edges  [u,  v )  of  sub(ei)  with  v  =  /oH'2(e,) 
(  see  Fig.  14  ). 

We  test  and  delete  repeatedly  as  described  above  until  we  find  an  edge  e,  that  does  not 
satisfy  AA  or  BB.  Then  we  can  construct  5,  recursively  from  sub^e,)  and  merge  it  into  G,_i. 
Since  no  edge  is  added  to  sub{e,)  during  the  construction  of  S,,  conditions  AA  and  BB  remain 
false  after  the  construction.  Thus,  the  resulting  graph  G,  will  be  planar,  and  no  /-normal  blocks 
will  interface. 

To  see  that  the  deleted  set  of  back  edges  is  minimal,  let  [u,  v]  be  an  edge  deleted  by  the 
above  algorithm,  and  add  it  back  to  G,.  If  [u,  v]  was  deleted  because  of  Condition  AA,  then 
low  1  (e,)  =  V  now,  and  Condition  AA  is  true  again.  If  [u,  v  ]  was  deleted  because  of  Condition 
BB,  then  /o>v2(e,)  =  v  now,  and  Condition  BB  is  true  again.  Notice  that,  in  the  latter  case,  the 
low  I  value  of  e;  has  remained  unchanged  since  the  deletion  of  [u,  v  ].  In  either  case,  G,  will  not 
be  /-planar. 

43.  Data  structures  and  running  time 

In  the  algorithm  described  above,  we  need  to  repeatedly  select  an  unprocessed  successor  of 
e  with  the  smallest  low  j  value,  and  the  low  i  values  of  tree  edges  are  constantly  changing.  There- 
fore we  maintain  a  heap  [13]  based  on  bw  i  values  of  the  unprocessed  successors  of  the  tree  edge 
e  currently  being  processed.  Since  the  algorithm  is  recursive,  we  actually  maintain  simultane- 
ously a  heap  of  unprocessed  successors  for  each  tree  edge  along  the  path  to  the  currently  active 
tree  edge.  The  total  size  of  all  such  heaps  is  0{m).  The  initialization  of  all  these  heaps  takes  a 
total  of  0{m)  time.  When  the  /owj  value  of  some  element  in  a  heap  increases,  we  modify  the 
heap  accordingly.  It  is  important  to  note  that  any  two  edges  in  active  heaps  are  unrelated;  thus 
deletion  of  a  single  attachment  can  modify  the  /owj  value  of  only  a  single  such  edge.  It  follows 
that  the  total  number  of  modifications  to  and  deletions  from  heaps  is  O  (m).  The  time  for  the  heap 
operations  is  0(}ogn)  time  per  operation,  for  a  total  of  0(,m\ogn)  time.  (Since  m  <  n^,  logm  = 
OQogn).) 

We  also  need  a  data  structure  for  the  back  edges  oi  sub(e)  so  that  the  following  operations 
can  be  done  efficiently: 

1.  delete  an  attachment  [u,  v  ]  of  e  with  v  =  low^(e)  or  v  =  /ow2(e); 

2.  maintain  the  low  \  and  low  2  values  of  e; 

3.  split  the  data  structure  into  several  pieces,  one  for  each  successor  of  e. 

One  easy  solution  that  meets  these  requirements  is  the  selection  tree  [8].  To  represent  a  set 
of  edges  £0  as  a  selection  tree  Tq,  we  store  edges  of  £0  inside  the  leaves  of  Tq  from  left  to  right 
in  increasing  order  (by  DFS  number)  of  their  tails.    Edges  with  the  same  tail  are  ordered 
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arbitrarily.  Each  internal  node  w  of  Tq  has  two  children  w.lchild  and  w.rchild.  Let  5h,  be  the  set 
of  edges  stored  in  the  leaves  of  the  subtree  rooted  at  w;  let  lb  =  min  {x:  [x,  y]€  S^]  and  rb  = 
max  { j: :  [x,  y ]  €  5„, ) ;  let  lowi  =  mm  {y:[x,  )/]&  S,^},  and  low2=  mm{{y:  [x,  y]e  S^  \ 
y^lowi  ]  u  (n-f-1}).  Then  the  four  values  lb,  rb,  lowi  and  I0W2  are  stored  in  the  four  fields 
w.lb,  w.rb,  w.low  \ ,  and  ■w.low2  of  w  respectively.  If  w  is  a  leaf  storing  the  edge  [x,  y],  then  w.lb 
=  x,  w.rb  =  x,  w./ow|  =  y,  and  W./0W2  =  n+\.  The  values  in  each  internal  node  can  be  computed 
from  the  values  in  the  children  ( in  constant  time  ). 

In  the  following  discussion,  we  will  refer  to  a  tree  by  its  root.  Let  ri  and  r2  be  two  selec- 
tion trees  representing  the  two  disjoint  sets  of  edges  E\  and  £2-  If  ui  <U2  for  !iil[ui,  vi]  e  £1 
and  [U2,  V2]  e  £2.  then  we  can  merge  ri  and  r2  to  get  the  selection  tree  for  £1  u  £2  in  0(1) 
time: 

procedure  merge(r^,  r^); 
begin    if  r,  =  nu// then 
return  r-^; 

end  if; 

if /-J  =  nu//  then 
return  r , ; 

end  if; 

r  :=  newnode  (); 

r.lchild  :=  r  1 ; 

r.rchild  :=  r2; 

r./i):=r,.Zi>; 

r.rfc  :=  r^-rb; 

r.towi  :=  mini/-]  ./otv,, /■2./oM'i ); 

r.low2  :=  min({ri  ./owi,  r^-low^,  r^  .low 2,  r^-low^]  -  [r.low^ )); 

return  r; 
end; 

Let  r  be  a  selection  tree  representing  a  set  of  edges  £o.  To  split  £o  into  two  sets  E\  = 
{  [m,  V  ]  e  £0  I  w  <  U;c }  and  £2  =  {  [m,  v  ]  e  £0  I  w  >  a^ } ,  we  split  r  with  respect  to  Ux  as  follows: 

procedure  split(r ,  u,); 
l)egin    if  tt,  <  r./ii  then 

return  [nulL  r]\ 
elseif    u,  >  r.rfc  then 

return  [r,  nult\\ 
else       [r;,  /-,]  :=  [r.lchild,  r.rchild]; 
if  «,  <  r,.rb  then 

['•n.'";2l:=.sp'"('-;.  ",); 
return  Ir,^, merge{r,2,  r,.)]; 
eise      [^,1,^,2]  ■.=  spla(r„  «J; 

return  [mcrge(r,,^,).'"r2]; 
end  if; 
end  if; 
end; 
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The  height  of  any  tree  that  results  from  splitting  a  tree  r  can  be  no  greater  than  the  height  of 
r.  To  select  and  delete  an  edge  {x,  v]  from  a  tree  r,  where  v  e  {r.lov^x,  r./ow2),  we  do  the  fol- 
lowing: 

procedure  deleie(r,  v); 
begin    if  r  is  a  leaf  then 

mark  the  back  edge  stored  in  r  as  'deleted'; 
return  null; 
else       (r,,  t,\  :=  [r.lchild,  r.rchild]; 

If  V  =  ri.low  1  or  V  =  r,.low 2  then 

return  merge{deleie(r,,  v),  r,); 
else       return  merge(r,,  delete(/-,,  v)); 
end  if; 
end  If; 
end; 

Assuming  the  input  graph  G  is  connected,  we  know  that  all  the  tree  edges  form  a  tree.  Let 
the  root  be  1.  For  technical  reasons,  we  add  a  dummy  edge  cq  =  [0,  1]  to  the  tree  edges.  To  get  a 
maximal  planar  subgraph  of  G,  we  just  construct  a  0-planar  subgraph  of  sub{eo)<  and  then  delete 
Co  from  it.  Initially,  we  construct  a  balanced  selection  tree  tree{eo)  to  store  all  the  back  edges  of 
G.  The  height  of  this  tree  is  0(logn).  The  lime  and  space  needed  to  initialize  treeico)  are  both 
0{m). 

When  we  begin  to  construct  a  maximal  /-planar  subgraph  for  a  tree  edge  e,  we  first  split 
tree(e)  into  several  pieces  tree(e\),  ...,  treeie^),  where  e^,  ....  e^  are  the  successors  of  e  not 
marked  as  'deleted'.  For  each  such  successor  e„  tree{ei)  is  a  selection  tree  representing  the  set  of 
back  edges  in  sub(ei),  and  can  be  obtained  as  follows.  If  e,  is  a  back  edge,  then  tree  (e,)  consists 
of  a  single  edge,  and  can  be  constructed  in  O (1)  time.  If  e,  is  a  tree  edge,  let  e,  =  [b,  c,],  and  let 
n,  be  the  number  of  descendants  of  c,.  It  is  well  known  that  a  back  edge  [u,  v  ]  is  a  descendant  of 
e,  iff  Ci<u  <  Ci  +  Hi.  Then  we  can  use  the  procedure  split  to  get  tree  (e,)  from  tree  {e)  in  O  (logn) 
time.  Since  there  are  at  most  n  tree  edges  in  G,  then  the  splitting  takes  takes  O(nlogn)  time  for 
the  whole  algorithm.  After  each  split,  the  total  size  of  the  trees  is  still  0(m). 

To  select  and  delete  an  attachment  [x,  v ]  of  e,,  where  v  e  {/owi(e,).  /oH'2(e,)},  we  execute 
delete  (tree  (ei),  v),  which  takes  O  (login))  time.  There  can  be  at  most  0(m)  such  invocations  of 
delete,  so  the  total  cost  for  executing  delete  is  0{m\ogn).  Given  the  selection  tree  tree(e,),  the 
values  (i)(e,),  lowiiei),  and  /ow2(e,)  can  be  computed  from  treeiei)  in  00)  time:  if  tree(e,)  is 
null,  we  just  set  these  values  to  n+\;  otherwise,  they  can  be  computed  from  treeiei).lowi  and 
tree(ei).bw2.  Thus,  the  total  cost  of  selection  tree  operations  is  O(mlogn). 

We  have  mentioned  that  the  total  cost  of  heap  operations  is  also  0(m\ogn).  The  other  costs 
of  the  algorithm  are  the  same  as  in  planarity-testing.  Thus  the  total  cost  of  our  maximal  planar 
subgraph  algorithm  is  O  (mlogn). 
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4.4.  The  complete  algorithm 

Now  we  summarize  our  maximal  planar  subgraph  algorithm.  We  take  a  connected 
undirected  graph  as  input,  and  convert  it  into  a  DFS  representation  G  =  (V,  T,  B).  At  the  same 
time,  we  compute  the  two  mappings  succ  and  N,  where,  for  each  e  e  T,  succ(e)  gives  the  succes- 
sor edges  of  e  in  increasing  order  of  their  tails,  and  for  each  v  e  V,  N{v)  gives  the  number  of  des- 
cendants of  V.  We  assume  that  there  is  a  dummy  edge  ^o  =  [0,  1]  such  that  succ(,eo)  gives  the  list 
of  tree  edges  leaving  the  root.  The  whole  preprocessing  takes  0(m)  time. 

We  summarize  the  maximal  /-planar  subgraph  algorithm  below. 

procedure  lplanar(e,  /); 
begin    let  e  =  [a,  b]; 

itee  fi  then 

return  [[[b],  (]]]; 

end  if; 

if  e  has  no  successors  then 
return  []; 

end  if; 

let  e , e^  be  the  successors  of  e  not  marked  as  'deleted'; 

split  tree  (e)  into  tree  (cj ) tree  (e^y, 

organize  Cj Cj  into  a  heap  based  on  their  /ow,  values,  with  the  smallest  one  on  the  top; 

let  e  I  be  the  edge  on  the  top  of  the  he^; 

delete  e  j  from  the  heap; 

1  alt{e):=lplanar{e^,l); 

delete  all  the  occurrences  of  b  from  the  lop  blocks  of  att  (e); 

y:=l; 

(•  :=  2; 

while  heap  is  not  empty  do 

let  e,  be  the  edge  on  the  top  of  the  heap; 

lf/ow,(e,)>ftthen 

delete  e,  from  the  heap; 

2  dummy  :=  Iplanar  (e„b); 
eiseif  Condition  AA  is  true  then 

V  :=  tree{ei).  low-^; 
while  V  =  tree(ei).  low^  do 

3  tree  (e, ) :  =  delete  {tree  (e, ),  v) ; 

end  while; 

if  e,  is  a  back  edge  then 

delete  e,  from  heap; 
eise  modify  heap; 
end  if; 
else  If  Condition  BB  is  true  then 

V  :=  /ree(e,).  low 2; 
while  V  =  tree(ei).  low  2  do 

4  tree(e;)  :=  delete(uee{e^).  v); 

end  while; 

else       delete  e,  from  the  heap; 

5  atl{ei)  :=  lplaruir(ei,  b)\ 
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6  merge  blocks  of  all  (e^)  into  one  block  B,; 

delete  all  the  occurrences  of  b  from  the  lop  blocks  of  ai((C{); 
If  I  =  2  then 

7  I)erfonn  Step  3  of  Alg.l; 

end  if; 

8  merge  fi,  into  all (e)  as  described  in  Step  4'; 

«:=i  +  l; 

end  if 
end  while; 
return  att(e); 
end; 

The  procedure  [planar  (e,  I)  implicitly  constnicls  a  maximal  /-planar  subgraph  of  sub  (e)  by 
deleting  a  minimal  set  of  back  edges.  The  parameter  /  specifies  where  the  blocks  of  afr(e)  are  not 
allowed  to  interlace  in  the  resulting  subgraph,  so  that  this  subgraph  can  be  used  to  build  a  larger 
planar  subgraph  without  further  deletion  when  we  process  the  predecessor  of  e.  For  the  initial 
call  where  e  =  eo.  we  have  /  =  0,  meaning  that  we  need  to  construct  a  maximal  planar  subgraph  of 
subieo).  In  the  recursive  calls  for  the  successors  of  e,  the  /  values  are  determined  as  follows. 
Since  no  /-normal  blocks  are  allowed  to  interlace  in  sub{e),  then  no  /-normal  blocks  are  allowed 
to  interlace  in  sub{e\)  either.  Thus  the  recursive  call  of  Iplanar  for  e\  (line  1)  has  the  same 
parameter  /  as  for  the  edge  e.  The  remaining  calls  for  ^2,  ....  Ck  Gine  5)  just  construct  maximal 
strongly  planar  subgraphs,  therefore  have  b  as  their  /  values.  Thus  when  we  merge  blocks  at  line 
6,  no  normal  blocks  of  att{e,)  interlace.  At  line  7,  we  merge  all  normal  blocks  of  att{e)  above 
low\{e2)  in  to  one  block;  at  line  8,  we  merge  B,  into  sub(e).  Because  of  the  deletions  at  lines  3 
and  4,  these  steps  can  be  performed  successfully  (without  any  further  deletion).  At  line  2,  sub(ei) 
is  detected  to  be  a  biconnected  component  and  is  processed  separately. 

To  compute  a  maximal  planar  subgraph,  we  simply  do  the  following: 

1.  Organize  B  into  a  selection  tree  treeieo); 

2.  Execute  lplanar{eo,  0); 

Then  T  <u  B  -  B'  gives  a  maximal  planar  subgraph  of  G,  where  B'  is  the  set  of  back  edges 
deleted  by  the  procedure  delete  in  the  preceding  algorithm. 

Remark  The  procedure  Iplanar  can  be  greatly  simplified  if  we  know  that  all  the  tree 
edges  of  G  are  on  a  same  cycle.  In  this  case,  the  low^  values  need  not  be  dynamically  main- 
tained: if  e  =  [a,  b]  is  a  back  edge,  then  low^ie)  =  b;  otherwise  lowyie)  =  1.  The  I0W2  values, 
which  are  used  only  for  testing  condition  BB  when  i  >  1.  need  not  be  maintained  either,  since  for 
/  >  1,  Ci  is  a  back  edge.  As  a  result,  the  selection  trees  are  no  longer  useful,  and  the  heaps  storing 
the  successors  of  tree  edges  can  be  replaced  by  lists  precomputed  as  in  Algorithm  1 .  With  these 
simplifications,  our  algorithm  gives  the  following  result  which  is  first  reported  in  [10]  by 
Kobayashi  et  at:.  A  maximal  planar  subgraph  can  be  constructed  in  linear  time  provided  that  a 
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Hamillonian  tour  of  the  graph  is  given.  ■ 

5.  Summary 

The  problem  of  drawing  graphs  in  the  plane  arises  naturally  in  circuit  layout.  Since  finding 
a  maximum  planar  subgraph  is  A'P-complete  [5],  a  maximal  planar  subgraph  seems  to  be  a  rea- 
sonable approximation.  Because  planarity-testing  can  be  done  in  linear  time,  it  is  easy  to  solve 
the  maximal  planar  subgraph  problem  in  O  (mn)  time:  start  with  a  graph  H  with  no  edge;  for  each 
edge  of  the  input  graph  G,  add  it  to  H  if  the  resulting  graph  is  planar,  and  reject  it  otherwise.  The 
resulting  graph  H  will  be  a  maximal  planar  subgraph  of  G.  However,  a  better  solution  seemed  to 
be  hard  to  find  for  a  long  time.  Jayakumar  et  al.  [9]  even  made  the  conjecture  that  "no  maximal 
planarization  algorithm  of  complexity  better  than  0(mn)  will  be  possible."  Our  O(m\ogn)  solu- 
tion disproves  this  conjecture,  as  does  the  method  of  Di  Battista  and  Tamassia  [3]. 

We  have  assumed  that  the  input  graph  to  our  algorithm  is  connected.  For  a  more  general 
graph,  we  can  find  a  maximal  planar  subgraph  by  applying  our  algorithm  to  each  of  its  connected 
components. 
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